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Rotating black holes in Einstein-Maxwell-Chern-Simons theory possess remarkable features, when 
the Chern-Simons coupling constant reaches a critical value. Representing single asymptotically flat 
black holes with horizons of spherical topology, they exhibit non-uniqueness. In particular, there 
even exist extremal and non-extremal black holes with the same sets of global charges. Both extremal 
and non-extremal black holes form sequences of radially excited solutions, that can be labeled by the 
node number of the magnetic gauge potential function. The extremal Reissner-Nordstrom solution 
is no longer always located on the boundary of the domain of existence of these black holes, and it 
neither remains the single extremal solution with vanishing angular momentum. Instead a whole 
sequence of rotating extremal J = 0 solutions is present, whose mass converges towards the mass of 
the Reissner-Nordstrom solution. These radially excited extremal solutions are all associated with 
the same near horizon solution. Moreover, there are near horizon solutions that are not realized as 
global solutions. 

PACS numbers: 04.40.Nr, 04.50.-h, 04.20.Jb 


I. INTRODUCTION 

In four spacetime dimensions vacuum black holes and black holes in the presence of an electomagnetic field possess 
very remarkable properties, as formulated by a number of important theorems. The black hole uniqueness theorem, 
for instance, states that asymptotically flat non-degenerate electrovac black holes are uniquely described by their 
global charges: the mass, the angular momentum and the electric (and magnetic) charge [1H4J. 

When going to higher dimensions, the Schwarzscliild solutions and the Reissner-Nordstrom (RN) solutions possess 
simple generalizations, while the Kerr solutions are generalized by the family of Myers-Perry (MP) solutions [§, Q. 
In D dimensions, these solutions possess N = [(D — 1)/2J independent spatial planes and thus also N independent 
angular momenta. In five dimensions both angular momenta are bounded for a given mass, a feature which does not 
hold for D > 5. 

The higher-dimensional generalizations of the Kerr-Newman solutions are not known in closed form. However, 
they have been studied perturbatively (Ml and numerically fl5l - [l8l ]. In odd dimensions, when all angular momenta 
possess equal magnitude, the angular coordinates factorize, leading to substantial simplifications of the perturbative 
and numerical studies. Moreover, for extremal black holes the near-horizon solutions have been investigated (l8l, Il9(. 
In the case of black hole solutions with equal magnitude angular momenta in odd dimensions, there are two branches 
of near-horizon solutions, the MP branch and the RN branch, intersecting at a critical point. Interestingly, for both 
branches, only a part of the near-horizon solutions also correspond to global solutions, thus not all Einstein-Maxwell 
(EM) near-horizon solutions have global counterparts. 

In odd dimensions, a Chern-Simons (CS) term can be added to the action, breaking the charge reversal symmetry 
for even N = (D — D/2. The resulting Einstein-Maxwell-Chern-Simons (EMCS) theories possess intriguing sets of 
black hole solutions [20--25]. In five dimensions for the special case of the CS coupling constant A = Asg, as obtained 
for minimal supergravity, the general set of charged rotating black hole solutions is known in closed form [22| . 

Carrying opposite angular momenta, the BMPV (20j black holes represent a subset of these solutions. Emerging 
from one of the charge symmetric extremal RN solutions, they form a branch of stationary extremal black hole 
solutions with vanishing horizon angular velocities. The magnitude of their (equal magnitude) angular momenta can 
be increased, while keeping their mass and charge fixed, until a singular solution with vanishing area is encountered. 
Recall, that in four dimensions stationary EM black holes with a non-rotating horizon are static. It is the presence 
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of the CS term, which allows for non-static black holes with vanishing horizon angular momenta. Here the frame 
dragging effects at the inside and outside of the horizon precisely cancel at the horizon in the case of the BMPV black 
holes [261] . 

When the CS coupling constant A is increased beyond its supergravity value, the set of black hole solutions acquires 
interesting new features. First, for Asg < A < 2Asg> counterrotation sets in, be., within a certain region of the domain 
of existence of the solutions, the horizon angular velocity and the angular momentum have opposite signs 24], 
As a consequence, extremal static black holes can become unstable with respect to rotation. Indeed, for fixed electric 
charge the mass can decrease with increasing magnitude of the an gula r momentum. Thus supersymmetry marks the 
borderline between stability and instability for EMCS black holes [23l . 126} . 

Next, for A > 2 Asg> the extremal RN solution remains part of the boundary of the domain of existence only for 
one sign of the charge. For the other sign of the charge, the extremal RN solution resides inside the boundary of the 
domain of existence. As boundary solution with vanishing angular momentum the static RN solution is now replaced 
by a set of two stationary solutions, whose global angular momenta vanish J23|. Thus EMCS theory allows for black 
holes with rotating horizon but vanishing angular momenta, where the contributions to the global angular momenta 
precisely balance. 

These two J = 0 rotating extremal solutions form only the lowest mass solutions of a whole sequence of excited 
J = 0 rotating extremal solutions, which possess an increasing number of radial nodes in one of the metric and one of 
the gauge field functions [25[ . For a fixed value of the charge, the mass of this sequence of solutions converges towards 
the mass of the corresponding extremal RN black hole. Since these excited extremal black hole solutions are located 
inside the domain of existence these solutions represent a new type of violation of uniqueness: there are extremal and 
non-extremal black holes with the same sets of global charges [25]. The violation of uniqueness among non-extremal 
black holes is of course also present [23}. 

When the full sets of extremal solutions are considered, an intricate web of branches arises. These branches of 
global black hole solutions can be compared with the branches of near-horizon solutions. The structure of the latter 
is much simpler. Comparison of the two sets shows, that there are near-horizon solutions, which correspond to (i) no 
global solutions, (ii) exactly one global solution, (iii) more that one global solution [25| . 

Here we give a detailed account of the astounding properties of EMCS black hole solutions in five dimensions, 
treating the CS coupling constant as a free parameter, a task which, to our knowledge, has not been yet considered 
in the literature. Special attention will be paid to extremal solutions, although some properties of non-extremal 
configurations will be discussed as well. The paper is organized as follows: in Section II we present the action, the 
Ansatze and the charges. We discuss the near-horizon solutions in Section III, where the near-horizon formalism must 
be employed with care because of the presence of the CS term. We present our numerical procedure for the global 
solutions in Section IV, discuss the boundary conditions and provide expansions for the functions. Our numerical 
results together with a comparison of global and near-horizon solutions are given in Section V. We end with a brief 
conclusion and outlook in Section VI. 


II. ACTION, ANSATZE AND CHARGES 

We focus our study on black holes in Einstein-Maxwell-Chern-Simons (EMCS) theory. Here we briefly review the 
action and the general set of equations of motion. We then present the appropriate Ansatze to obtain rotating black 
holes with equal magnitude angular momenta. Next we recall the general formulae for the global charges and the 
horizon properties of these black holes, as well as their scaling symmetry. 


A. Einstein-Maxwell-Chern-Simons action 


The action of Einstein-Maxwell-Chern-Simons theory in five dimensions reads 


1 = 


1 


167tG,5 


d 5 a 




12^3 


E ^fh AtiKaFfii 


( 1 ) 


where R is the curvature scalar, G 5 is Newton’s constant in five dimensions, A M is the gauge potential with field 
strength tensor F^ u = — d^A^, and A is the Chern-Simons (CS) coupling constant. In the following we employ 

units such that 167 tG 5 = 1 . Note, that for A = Asg = 1 the action corresponds to the bosonic sector of minimal 
supergravity. For A = Aem = 0 the action corresponds to Einstein-Maxwell theory. 

Variation of the action with respect to the metric leads to the Einstein equations 




( 2 ) 
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where the stress-energy tensor is given by 


T„ v = F w F v p - \g^F pG FP°. (3) 

Variation with respect to the gauge potential leads to the Maxwell-Chern-Simons equations 

V v F^ u + -^=e^ a ^F va F Pl = 0. (4) 


B. Ansatze 

We consider stationary black holes, which represent generalizations of the five-dimensional Myers-Perry solutions 
P] to EMCS theory. Hence these black holes possess a spherical horizon topology. Their Killing vectors are 

£ = dt, V(i)=d vi , rj (2) = d V2 , (5) 

with t the time coordinate and i pi, 2 angular directions. With each azimuthal symmetry an angular momentum J(fc), 
k = 1,2, is associated. In general, both angular momenta are independent. 

Here we restrict to the case where both angular momenta have equal magnitude, |«/(i)| = | J( 2 )| = \J\. The spacetime 
then becomes a cohomogeneity-1 manifold, with an enhancement of the isometry group from R t x U(l) 2 to Rt x U{ 2). 
Hence the angular dependence of the metric and the gauge potential can be explicitly given. Then the metric may 
be parametrized by the Ansatz 

ds 2 = —Fo(r)dt 2 + F\{r)dr 2 + F 2 {r)d9 2 + F 3 (r) sin 2 9 ( Eid<pi — W(r)dt ) 2 

+F 3 (r) cos 2 9 (e 2 dip 2 - W{r)dt) 2 + (F 2 (r) - F 3 (r)) sin 2 0cos 2 6(eidipi - e 2 dip 2 ) 2 , ( 6 ) 

where 9 £ [0, 7 r/ 2 ], ipi £ [0, 2n\ and ip 2 £ [ 0 , 27 t], and e*, = ±1 denotes the sense of rotation in the fc-th orthogonal 
plane of rotation. Note also that the line element Eq. f 6 ]) still possesses a residual metric gauge freedom. 

The corresponding Ansatz for the gauge potential is given by 

A M dx M = ao(r)dt + a v (r)( sin 2 9e\dyi + cos 2 9e 2 dtp 2 ). (7) 


C. Known solutions and parametrization used in the numerics 


For the general family of EMCS solutions within the above Ansatze, only the following special cases are known in 
closed form: 

(z) Vacuum rotating black holes (Myers-Perry solution). This solution is valid for every value of A, since the Maxwell 
and the Chern-Sinrons terms vanish identically. The functions which enter the line element Eq. ([ 6 ]) are 


F o(r) 

F 3 {r) 


1 _ M _i_ 3 J 2 

1 6tt 2 t- ' 8 7i - 2 Mr 4 

1 4- 3J2 

x ' &7r‘ 2 Mr 4 


Fi(r) = 


^ -I- 3J2 

^ ' 8-rr' 2 Mr 4 


F 2 (r) = r 2 , 


3 J 2 


87 r 2 Mr 4 


), W(r) = 


J 


4 7r 2 r 4 (l 


3 J 2 ' 

8ir 2 Mr 4 > 


( 8 ) 


while a v (r) = a 3 (r ) = 0. Also, M and J are the mass and the angular momentum of the solutions. 

(ii) Charged static black holes (Reissner-Nordstrom solution). Again this solution is valid for every value of A. In 
the static case, the magnetic component of the gauge field vanishes, a v (r) = 0, and the Chern-Simons term does not 
contribute. The expressions of the other functions which enter the solution are 


*b(r) 


1 _ M Q 2 

Fi(r) 67 r 2 r 2 487 r 2 r 4 ’ 


F 2 (r) = F 3 {r) = r 2 , 


W{r) = 0, and ao(r) 


Q 

47 r 2 r 2 ’ 


(9) 


with Q the electric charge. 
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(in) The general set of charged rotating solutions for A = Asg = 1 [23] with the extremal BMPV black holes 
as a special case. These black holes have 


l _ 2 (p-q) , 2 j 2 p+q 2 

Fo(r) = — ^=^-=5^, F\(r) = 


1 + 


1 


1 _ 2(p-q) | 2j 2 p+t/ 2 
1 ^2 I ^4 


F 2 (r) = r 2 , 


( 10 ) 


r, ^ , 2 j 2 p j 2 ? 2 \ 2 wt \ 2 P~q~h j , ^ ^ V3q 

p 3 (r) = (i + - ^-)r , W(r)= i + ^^ j2q2 and a*(r) =-a 0 (r) = —, 


with „ - M+V3Q ■ _ .3 J _ Q 

witn p - r2n2 , j - 2M+VSQ , q - 4V5t2 ■ 

In the numerics, we have found it useful to use quasi-isotropic coordinates by choosing a metric gauge with T^r) = 
F\(r)r 2 and take F 0 (r) = f(r), Ti(r) = jjfj, p 3 (r) = y|yyr 2 , W(r) = ^r-. This results in the following expression 
of the line element in terms of four unknown functions 


m(r) 


n(r) 


ds 2 = — f(r)dt 2 + (dr 2 + r 2 d 0 2 ) + —(r sin 2 6 Sidipi — 


f(r) 

originally proposed in 


f(r 

n(r) 9 o „ ( , w(r) 

+ r 2 cos 2 6 £ 2 dip 2 - 


dt 


f(r) 

m(r) — n(r) 


w(r) 


dt 


f ( r ) 


r 2 sin 2 9 cos 2 9(sidipi — £ 2 ^ 2 ) , 


(ID 


D. Global charges 

We here consider stationary asymptotically flat black holes. Their total mass M and angular momentum can 
be obtained directly from the Komar expressions associated with the Killing vector fields 


M = 


’si. 


( 12 ) 


J(k) — [ P(k) ) (13) 

JS 3, 

where a^ lM2M3 = e /ilAt2At3 p Cr V p ^ cr , and ^(fe)p lM2M3 = £ / j. 1 p, 2 p 3 p^ Pr l(k)- For equal-magnitude angular momenta J (fc ) = 

Skd , ^ — 1 , 2 . 

The electric charge Q is given by 


< i4) 

where F /il/J2 p 3 = e MlM2 M 3 pa F pcJ . The magnetic moment /x mag is determined from the asymptotic expansion of the 
gauge potential a v (r), see Section IV B. The gyromagnetic ratio g is then obtained from the magnetic moment p, mag 
via 


QJ 

Mmag ~ 9 2M ' 


(15) 


E. Horizon Properties 

The black hole horizon is located at r = rn and rotates with angular velocity Hr ■ It is a Killing horizon, where the 
Killing vector £ = d t + fln(£id Vl + £ 28 ^ 2 ) becomes null and orthogonal to the other Killing vectors 


(C 2 )lw = 0> (C‘dt)|«=0, (C • d<pi)\-H = 0, (C ’ <9<p 2 )|m = 0. 


( 16 ) 


Note, that these expressions impose conditions on the metric functions. 
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In particular, for the parametrization Eq. m employed in the numerics, the metric function f(r) must vanish at 
the horizon. 


f(ra) = 0. 


while the horizon angular velocity is given by 


Oh 


'Ch 


(17) 


(18) 


We focus our study on extremal black holes. In quasi-isotropic coordinates, the event horizon of extremal black 
holes is given by ru = 0. Note that in this case, the horizon angular velocity is obtained from 


Oh 


lim 


^(t-h) 


o/(r)|«. 


The area of the horizon Ah is given by 


A n= f \/lff (3) l = r H A (S 3 ) lim 
Jh v r ^ ri 



(19) 


( 20 ) 


and the surface gravity k reads 

« 2 = ~(VC) 2 |«= lim t 

z r->r H (r — rH)v m 

For extremal black holes the surface gravity vanishes, k = 0. 

The horizon mass Mh and horizon angular momenta Jn(fe) ar e given by 

M h = [ a, 

Z J-H 

Jn(k) = / P{k) ■ 

Ju 

For equal-magnitude angular momenta Jn(fc) = £kJn, k = 1,2. 

The electric charge Q can also be determined at the horizon, 

The horizon electrostatic potential <I>h is defined by 

HIh = C^A^H = (^0 d - • 

Note that Th is constant at the horizon. 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


F. Scaling Symmetry 


The EMCS solutions have the following scaling symmetry [23] ( with r > 0 an arbitrary parameter) 

M = t 2 M, Ji = T 3 Ji, f H = rr H , 0 H = O h /t, R = k/t, (26) 

etc. Let us therefore introduce scaled quantities, where we scale with respect to appropriate powers of the mass. 
These scaled quantities include the scaled angular momentum j = J/M 3 / 2 , the scaled charge q = Q/M , the scaled 
area an = Ah/M 3 / 2 , the scaled surface gravity R = kM 1 / 2 , and the scaled horizon angular velocity Oh = OhM 1 / 2 . 
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G. Smarr Formula 


The Smarr mass formula for EMCS black holes with two equal-magnitude angular momenta reads 


M = 3kAh + 3U H J + j 


or in terms of scaled quantities 


1 = 3koh + 3U H j + ■ 

Note that in the extremal case, the surface gravity vanishes. Moreover, in the extremal case 

Mu = Sfindn. 

EMCS black holes satisfy the first law [26j 


dM = 2 k cL4jj + 2fljj dJ + dQ. 


(27) 


(28) 


(29) 

(30) 


III. EXTREMAL CONFIGURATIONS: NEAR-HORIZON SOLUTIONS 

Before analyzing the numerical solutions, some analytical understanding of the properties of extremal EMCS black 
holes can be achieved by deriving the near-horizon solutions in the entropy function formalism f27l - i-29| . Employing 
this formalism we obtain semi-analytic expressions for the entropy as a function of the electric charge and the angular 
momentum. We note, however, that the near horizon formalism needs special care, when it is employed in the presence 
of a CS term, as discussed extensively in the literature (see e.g. il). 

In the following we first review the general entropy function formalism in the presence of a CS term. Subsequently, 
we apply the formalism to construct the near-horizon solutions of the extremal black holes in (i) pure Einstein-Maxwell 
theory, (ii) the bosonic sector of D = 5 supergravity, and (iii) EMCS theory with general CS coupling constant A. 


A. Entropy function formalism 


To apply the entropy function formalism to obtain the near-horizon geometry of extremal EMCS solutions, we make 
use of the Ansatz 


d r 2 

ds 2 = vi(— - r 2 dt 2 ) + V2[4:d9 2 + sin 2 29(£2dtp2 ^ £idipi) 2 ] 

+ V 2 il[eidipi + £ 2 difi 2 + cos26»(£ 2 d^ 2 — £idip±) — ardt] 2 . 

Note the shift of the radial coordinate r —> r — rn- Thus the horizon is located at r = 0. 
For the gauge potential we employ the Ansatz 

A = —( p+pa)rdt + 2 p (sin 2 9£\dip\ + cos 2 ^cfy^)- 


(31) 


(32) 


Here the parameters V 2 , n, p and p are constants, which satisfy a set of algebraic relations, that follow within 
the near-horizon formalism |27 j29|. 

In the near-horizon formalism, the entropy is obtained from the extremum of the entropy function 


S = 2tt( 2 aJ + pq — h), 

where h is an action functional, that depends on the constants of the near horizon Ansatz 


(33) 


h(a,vi,V 2 ,r},p, p) = J dddipid^y/^C 


8tt 2 


:(9'u|t 7 3 / 2 a 2 — 9v2vfif/ 2 — 36v 2 y/rjvi + 36vf^/rjv2 + 36v%yA]p 2 


9V\^/V2 

72v 2 sfrjppa + 36v 2 ^/rjp 2 a 2 — 36vly/rjp 2 — 32\V3p 2 pv\y/v 2 — 32\V3p 2 ‘avi\/v 2 


(34) 


and q is related to the electric charge (see Eq. (T43l) b 
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In the usual analysis, one proceeds by taking the derivative of the functional h with respect to the constants 
introduced in the Ansatz. These then yield the set of equations to be solved and the conserved charges J and Q. 
However, in the presence of a Chern-Simons term in the action the analysis needs to be modified [30]. 

To obtain the near-horizon geometry we can directly solve the Einstein equations and the Maxwell equations. The 
set of Einstein equations is equivalent to taking 


dh 

dv\ 


= 0, 


dh 
dv 5 


= 0, 


dh 

dr] 


= 0, 


(35) 


i.e., these equations continue to hold also for EMCS theory. Variations of the functional with respect to p, however, 
do not take into account the Chern-Simons contribution. Therefore the Maxwell equations must be used when A ^ 0. 
Simply setting dh/dp = 0 would lead to a wrong result in this case. 

The Einstein equations then yield the following algebraic relations 


v 2 = Vi, 


VVl = ~3 


4 (p — p + pa) (p + p + pa) 


a 2 — 1 


Vi = 


2 a 4 p 2 — p 2 + 2 a 3 pp — 4 ppa + a 2 p 2 — 2 p 2 


3 a 2 -l 

and the Maxwell equations lead to 

3av^ 2 ^/rjp + 3a 2 v^ 2 ^/rjp — AXVSpviV 2 p — 4:XV3p 2 ViV 2 a — 3pvly/v 2 y/r] = 0 . 


(36) 


(37) 


Note, that we now have 4 algebraic relations for 6 unknown constants. Thus we have two independent parameters. 
For those we may choose the angular momentum J and the electric charge Q. 

In the usual near-horizon formalism, the angular momentum J is obtained by taking the derivative of the action 
functional h. with respect to the associated constant a, and the electric charge Q is obtained by taking the derivative 
with respect to p. To see where modifications arise because of the CS term we now recall the expressions for the 
angular momentum and the electric charge as Noether charges [dpi . 1 3214341 . 

The Noether charges can be obtained by integrating the corresponding charge densities [30] 


or=- 


(V a £ M - W£ a ) - 4 (CA r ){V=gF a » + 


2X 

37I f 


^"PApF^) 


over the 5' 3 -sphere. When taking the Killing vector or ? 7 ( 2 ), the angular momentum J is obtained 

/» 3/2 5/2 

J= / dOdipid^Q*/ = 647r 2 —— y/rjp(p + pa) + 16ir 2 — — rf^ 2 a - —\f^Ti 2 p 3 X. 

J ^ v\ vi 9 


Since the very same result is obtained from the equation 

dh 

da 


= J, 


(38) 


(39) 


(40) 


the usual near-horizon relation remains valid for the angular momentum. 

However, this is not true for the electric charge Q. For Q we need to consider the charge density 


QZ = 4 


V~gF c 


x 


^ p ApF ; 


vp 


When we integrate this charge density over the S^-sphere we obtain the charge Q 

- 3/2 128tt 2 v / 3 


Q = / d8dipidip 2 QM = -647T 


tr - R/ '" 2 ^ 2 —^/ri(p+pa)+ vu Xp 2 . 

f V\ 


(41) 


(42) 


1287t 2 v / 3 2 

-Ap 2 , 


Note, that this is equivalent to 


dh 

dp 


q = -Q 


9 


(43) 
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and thus dh/dp ^ —Q, in contrast to a naive application of the entropy function formalism. Consequently, the 
extremization of the entropy functional must be appropriately modified to obtain the proper set of equations and 
charges. 

We are also interested in obtaining expressions for the horizon properties. The horizon angular momentum can be 
calculated from the standard Komar formula 

r 5/2 

J H = / - V’V) = 167T 2 ^_77 3 / 2 a. (44) 

J Vi 

Finally, for the horizon area we find the following expression 

Ah = J d9d(f>idtp2\l | det(</( 3 ))| = lbTr 2 ?;^ 2 ^. (45) 

Unfortunately, for a generic nonzero CS coupling constant A, A ^ Asg, it is not possible to give an explicit expression 
for S = ^4r/ 4G5 as a function of Q and J. A straightforward numerical analysis of the algebraic relations reveals a 
rather complicated picture, with several branches of solutions. We are now going to describe several cases in detail. 


B. Einstein-Maxwell: Aem = 0 

In the pure EM case the CS term vanishes, since A = 0, and the above relations reduce to those obtained in the 
usual entropy formalism. The EM case is interesting, since it allows for two branches of near-horizon solutions fT/l-flot . 
The near-horizon branch starting from the Myers-Perry solution is given by 


vi = v 2 , 

2(i'2 - 2 p 2 ) 

1 =-, 

V2 

P = 0, 

a = — 1, (46) 

J = 327r 2 U2 /2 v / 2-4p 2 /t>2, 

Q = 647T 2 p v / ^2 a/ 2 - 4 p 2 /v2, 

where J and Q are the angular momentum and the electric charge, respectively. The horizon angular momentum Jh 
and horizon area Ah are given by 

Jh = 167r 2 u2 /2 (2 - 4p 2 /u2) 3 ^ 2 , 

A h = 16tt 2 vI /2 v / 2 - 4 p 2 /v 2 = J/2. (47) 

Hence, along the full Myers-Perry branch, the entropy is always proportional to the angular momentum and indepen¬ 
dent of the charge. 

A second near-horizon branch starts from the Reissner-Nordstrom solution, and is given by 


Vi = v 2 = 


4(a 2 + l)p 2 
3(a 2 - l) 2 


V = 


P = 


a 2 + r 

pa 


* 2 - r 


512v / 37T 2 y0 3 Q! 

= 9(a 2 - l) 3 
128v / 3vr 2 /o 2 
Q = 3(a 2 - l) 2 ' 


(48) 
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It has horizon angular momentum and horizon area 

128^71 - 2 p 3 a J 

~ 9(a 2 - l) 3 “ 4 ’ 

A H = 128 ^ 2 / ( 1 Q ,3 +1) = J 2 Q~ 3/2 + ^V 2 - (49) 

9(ar — l) 0 ’ 2 487T 

Along this branch, the area is not proportional to the angular momentum. Instead, the horizon angular momentum is 
proportional to the total angular momentum. Interestingly, only parts of the near-horizon branches are also realized 
as global solutions [SHI- 


C. D = 5 supergravity: Asg = 1 

In the case of D = 5 supergravity, Asg = 1, the global solutions are known j22j , with the BMVP solution representing 
the special case of supersymmetric and thus ergofree solutions [2Cj ■ It is straightforward to obtain the near-horizon 
solutions in this case. They can be expressed as [30| 


vi = v 2 = fi/4:, 

_ , _ f_ 

71 1 fi 3 ’ 


a = 


P = 


vV - j 2 ’ 
VSfJ 2 

W m 3 - j 2 ’ 


V3j 

p= --» 

J = —An 2 j, 

Q = 8 \/Stt 2 fx, 


(50) 


where the parameters /r and j essentially describe the electric charge and the total angular momentum. 
The horizon angular momentum and the horizon area are given by 


r 2 K 2 j , 3 , 2 . 

J H = —s“(M - J ). 

fj, 

Ah = 2tt 2 v / |m 3 - J 2 | 


V2 

487T 


V3Q 3 - 2887r 2 J 2 


(51) 


This shows, that there are again two branches of near-horizon solutions. The first branch has J 2 >- ±=-Q 3 while 

3v37t 

the second branch has J 2 < — —4= Q 3 . The second branch thus corresponds to the BMPV branch. In terms of the 

3 v 37 t 

global solutions, this branch is ergo-region free with vanishing horizon angular velocity. In contrast, the first branch 
corresponds to the ordinary branch which possesses an ergo-region. At the matching point of both branches the area 
vanishes. 


D. Generic values of A 

Let us now consider the near-horizon solutions for generic values of the CS coupling constant A. Note, that the 
solutions cannot be given explicitly, except for the previous two cases. Also, for A 4 0 the CS term breaks the charge 
reversal symmetry Q —¥ — Q for spinning solutions. Therefore, we need to consider positive and negative Q separately 
for finite values of A. 

Let us first consider near-horizon solutions with positive Q. When A is increased from zero, the branch structure of 
the pure EM case is lost. Instead a new branch structure emerges, which is generic for 0 < A < 1.91. It resembles the 
supergravity case, which is included as a particular case. After a small transition region 1.91 < A < 2 a new generic 
branch structure arises for A > 2. 
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(a) (b) 


FIG. 1: Near-horizon solutions: The horizon area Ah (a) and the horizon angular momentum Jh (b) versus the angular 
momentum J for CS coupling constant A = 0, 0.5, 1, 1.5, 2, and 3 and charge Q = 4. The symmetric solutions obtained for 
J —> — J are suppressed. 


Let us now demonstrate this A-dependence in detail. In Fig. |1 (a)] we exhibit the horizon area An versus the angular 
momentum J for several values of A, where A is increased from zero to three and the charge is fixed at the positive 
value Q = 4. Note, that the set of symmetric solutions obtained for J —> — J is not exhibited here. It is neither 
exhibited in Fig. |l(b)| which shows the horizon angular momentum Jh versus J for the same set of parameters. 

Keeping in mind that a set of symmetric solutions is obtained for J — J, we focus for the moment the discussion 
on the solutions exhibited in the figure. For A ^ 0 and positive Q , there exists always a solution with vanishing area. 
The angular momentum J of this solution is finite. For fixed Q its angular momentum decreases monotonically with 
increasing A, as seen in Fig. 1(a) Fig. |l(bj| demonstrates, that the horizon angular momentum changes smoothly at 
these solutions with vanishing area. 

In the interval 0 < A < 1.91 there are two near-horizon branches. The small-J branch extends from the extremal 
Reissner-Nordstrom solution at J = 0 to the solution with vanishing area, which has the maximal value of J along 
this branch. The large-J branch, on the other hand, extends from the solution with vanishing area and the minimal 
value of J along this branch to solutions with arbitrarily large values of J (Myers-Perry limit). Thus the structure of 
the near-horizon branches is essentially the same as the one found for the supersymmetric value Asg = 1 , except that 
for Asg = 1 the small-J branch is ergo-region free. 

Let us next address the transition region 1.91 < A < 2. Here the branch structure changes. This is illustrated in 
Figs. |2(a)| and [2(b)] Note, that the set of symmetric solutions obtained for J —> —J is again not exhibited here to 
have more clarity in the figures. 

When the horizon area An is considered as a function of the angular momentum J, 
swallowtail bifurcation is seen for any value of A in the transition region 1.91 < A < 


as shown in Fig. 2(a) 
Hence 


2. Hence, there are solutions 
with the same value of J but different values of the horizon area or the horizon angular momentum in the transition 
region. The reason for the emergence of the swallowtail bifurcation becomes clear by inspecting Fig. |2(b)| where the 
angular momentum is shown versus the near-horizon parameter v\. Here an inflection point arises at A = 1.91, and 
entails the occurrence of a local minimum and a local maximum for larger values of A. 

With increasing A the endpoints of the two cusps of the swallowtail bifurcation move to smaller values of J and 
separate from each other. This has two interesting effects. The first effect is that, at the particular value of A = 1.9425, 
the upper cusp endpoint reaches the value J = 0. The associated black hole solution thus corresponds to an extremal 
stationary black hole, that is not static, i.e., it differs from the extremal static Reissner-Nordstrom solution, that is 
also present. (Of course, by the symmetry J —> — J there are two such stationary nonstatic solutions. These have the 
same area but differ in the signs of their horizon angular velocities and horizon angular momenta.) 

As seen in Fig. |2(a)| for 1.9425 < A < 2 the upper cusp endpoint is at negative values of J and thus the number of 
nonstatic J = 0 solutions is doubled. With increasing A, the upper cusp endpoint continues to move to lower values 
of J, and so does the lower cusp endpoint. However, for any value of A, the extremal Reissner-Nordstrom solution 
also resides at J = 0, and is thus approached by the lower cusp endpoint. This leads to the second interesting effect, 
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(a) (b) 

FIG. 2: Near-horizon solutions: (a) The area Ah versus the angular momentum J for values of the CS coupling covering the 
vicinity of the transition region near the critical value A cr = 2. The value of the charge is Q = 8%/37t 2 . The symmetric solutions 
for J —» — J are not shown, (b) The angular momentum J versus the parameter vi for the same sets of solutions. 


namely the disappearance of the lower cusp at the critical value A cr = 2, and thus the disappearance of one of these 
(two pairs of ) J = 0 solutions. _ 

In particular, as seen in Fig. |2(b)| when the critical value A cr = 2 is approached, the local minimum in J deepens 
further, while the local maximum merges with the endpoint corresponding to the static solution. From the point of 
view of the full set of solutions, i.e., including the symmetric J —> — J set of solutions, at the critical value A cr = 2 
two swallowtail bifurcations merge to form a single swallowtail bifurcation. Thus at A cr = 2 in total only one pair of 
nonstatic J = 0 solutions remain together with the static Reissner-Nordstrom solution. 

id |3(b)| 


Beyond the critical value A cr = 2 we observe a new generic branch structure, exhibited in Figs. 3(a) 


anc 


These 

figures show the full set of near-horizon solutions, including the symmetric set of solutions obtained for J —> — J as 
well as the solutions for negative values of the charge. In the figures the value A = 5 is chosen for the CS coupling 
constant. 



(a) 


(b) 


FIG. 3: Near-horizon solutions: The horizon area Ah (a) and horizon angular momentum Jh (b) versus the angular momentum 
J for CS coupling constant A = 5. The value of the charge is \Q\ = 1. 


For positive Q these near-horizon solutions possess three distinct branches, while they feature only a single branch 
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for negative Q. Because of the symmetry J —> — J, all these branches have mirror branches in Fig. 3(a) Considering 


first positive Q , the small-J branch starts from the extremal Reissner-Nordstrom solution and extends until the cusp 
Cnh is reached. The intermediate-J branch extends from the cusp C n h to the solution with vanishing area. From 
there the large-J branch extends towards solutions with arbitrarily large values of the angular momentum, reaching 
in the limit J —> oo the Myers-Perry solution. 

Clearly, beside the static Reissner-Nordstom solution there are two symmetric J = 0 solution, which are nonstatic. 
As seen in Fig. |3(b)| their horizon angular momenta are finite and opposite to each other. These nonstatic J = 0 
solutions are only present in the set of positive Q solutions. Likewise the solutions with vanishing area are only present 
in the set of positive Q solutions. In contrast, the negative Q solutions show a rather unspectacular behavior. Their 
single branch connects directly from the Reissner-Nordstrom solution to the Myers-Perry solution for any finite value 
of the CS coupling constant A. 

Concluding, we have seen that the branch structure of the near-horizon solutions depends significantly on the value 
of A and on the sign of the charge Q. In the next section we discuss the numerically obtained sets of global solutions 
and compare with the corresponding sets of near-horizon solutions. We show that for negative Q, all near-horizon 
solutions correspond to global solutions. In contrast, when positive values of Q are considered, not all of the near¬ 
horizon solutions are realized globally. Moreover, other near-horizon solutions may correspond to more than a single 
global solution. 


IV. GLOBAL SOLUTIONS: NUMERICAL PROCEDURE 

Having obtained analytical expressions for the extremal solutions in the near-horizon formalism, we are now inter¬ 
ested in obtaining global solutions, both extremal and non-extremal. As discussed in Section II, for generic values of 
the Chern-Sinrons coupling constant A and of the global charges, numerical integration of the differential equations 
seems necessary. 

In the following we discuss the numerical procedure, the boundary conditions and the expansions of the functions 
employed to obtain generic EMCS black holes in five dimensions. 


A. Differential equations and boundary conditions 


Considering variation of the Lagrangian with respect to the ag(r) component of the gauge field we can obtain a 
first integral of the system in terms of the electric charge of the black hole. 



w(r) d 
r dr 


a v (r) + 


Hr ) 3 / 2 

yj m(r)n(r) 


4 

3 


v / 3Aa l/ ,(r) 2 



(52) 


This differential equation is compatible with the Einstein and Maxwell equations. Combining this first-order differ¬ 
ential equation with the system of differential equations, we obtain a minimal system of differential equations for the 
EMCS black hole solutions: four second-order differential equations for /(r), m(r), uj(r ), and a v (r), one first-order 
differential equation for ro(r), and one first-order differential equation for ao (r), which is decoupled from the other 
differential equations. 

For the numerical calculations, we have found it useful to introduce a compactified radial coordinate. For the 
non-extremal solutions we take the compactified coordinate to be x = 1 — r^/r. In the extremal case we employ 
x = jzrp- (Note, that we are using an isotropic radial coordinate r, so rn = 0 in the extremal case.) We employ 
a collocation method for boundary-value ordinary differential equations, equipped with an adaptive mesh selection 
procedure jUJ. Typical mesh sizes include 10 3 — 10 4 points. The solutions have a relative accuracy of 10 _1 °. The 
estimates of the relative errors of the global charges and the magnetic moment are of order 10~ 6 , giving rise to an 
estimate of the relative error of the gyromagnetic ratio g of order 10 -5 . 

To obtain asymptotically flat solutions the metric functions should satisfy the boundary conditions at infinity 


f\r=oo — ITllr-oo — TL |r—oo — 1 : ^|r=oo — 0 • 


(53) 


For the gauge potential we choose a gauge in which it vanishes at infinity 


U()|r=oo — &Lp\r=oo — 0 ■ 


Requiring the horizon to be regular, the metric functions must satisfy the boundary conditions 


(54) 


f\r=r H — ?Tl\r=rii — Tl\ r —rn — ^ 5 ^|r=rn — 7‘H^H ^ 


( 55 ) 
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where Hh is the horizon angular velocity, Eq. m- The gauge potential satisfies at the horizon the conditions (l25l) 


C M Au| r=rH = $ h = (ao + ^Ha v )| r=m 
with the constant horizon electrostatic potential $h- 


dap 

dr 


r—r h 


= o, 


(56) 


B. Expansions 


From the asymptotic expansion of the metric functions and the gauge field functions we can extract various param¬ 
eters of the black hole: 


/-> 1 - - 
J f 

J 


M 


4 7r 2 r 3 


a o 


m —> 1 — 

4t r 2 r 2 ’ 


M 


, n —> 1 — 

M mag 


M 


I2ir 2 r 2 


4w 2 r 


(57) 


where M is the mass, J the angular momentum and /j, ma g the magnetic moment. From these charges we can calculate 
the black hole gyromagnetic ratio g = 


2M/ln 


qj ; 

Next we consider the expansion of the functions near the horizon. First we present the expansion for non-extremal 
black holes. Once the previous conditions are imposed, it can be seen that: 

f(r) = / 2 (r-r H ) 2 +o((r-r H ) 3 ), 
m(r) = m 2 (r-r H ) 2 +o((r-r H ) 3 ), 

n(r) = h(r — r H ) 2 + o((i r H ) 3 ), (58) 

u>(r) = Hh + o(r — r H ), 
a 0 (r) = a 0 , 0 + o((r - r H ) 2 ), 
a v (r ) = ap fi + o((r - r H ) 2 ). 

These constants are implicitly related to three global charges (M, Q and J) in the non-extremal case. 

In the case of extremal black holes, the radial dependence of the functions in the vicinity of the horizon changes as 
follows: 


f{r) 

m{r) 

n(r) 

w(r) 

ao(r) 

M r ) 


Ur 4 + / Q r(“ +4 > + o(r 6 ), 
m 2 r 2 + mpr(P +2 ' 1 + o(r 4 ), 
hr 2 + l 1 r ( ' 1+2 ' > + o(r 4 ), 

+ w 2 r 2 + o(r 3 ), 

Oo.o + ao,\r S + o(r 2 ), 
a v ,o + lk, M r M + o(r 2 ). 


(59) 


Since we are considering the extremal case, all these constants are implicitly related to two global charges (Q and J 
for example). The coefficients a, /?, 7 , A and g can be non-integer. In fact we have 


0 < a < 2 , 0 < /3 < 2 , 0 < 7 < 2 , 

0 < <5 < 2 , 0 < fj, < 2 , 0 < is <2 . 


(60) 


In the pure Einstein-Maxwell case (A = 0) the expansion contains similarly non-integer exponents. This feature is 
found in both EM branches (be., in the MP branch and in the RN branch). 

Thus the expansion at the horizon in general leads to non-integer exponents. We observe that these exponents (a, 
/?, 7 , A and g) can become lower than one. This could lead to a divergence of the first derivative of a reparametrized 
function. Hence we have to be very careful with the reparametrization we use for the numerics. We reparametrize 
the functions in the following form: 

m = m(r) 
fi = n{r) 

f = f( r )/x 2 
Cj = (1 — x)~ 2 uj(r) 
dp = x 2 ap(r). 


( 61 ) 
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We then integrate the functions m,n,f,ui and a v . The advantage of using these reparametrized functions is the 
following. Since the exponents (a, /?, 7 , A and g) are always greater than 0, and the order of the differential equations 
is lower or equal than 2 , there is no problem with diverging functions or derivatives in the numerics. 

Note, that all the redefined functions except for Cj now start with an a; 2 -term in the compactified coordinate 
x = r/(r + 1). (The reparametrization of u> is not related to the expansion at the horizon. It is done in order to be 
able to fix the angular momentum by a boundary condition). 

The boundary conditions must guarantee that the functions satisfy the horizon expansion ( m,n,f and a^, start 
as x 2 ). At the horizon, x = 0, we can impose the angular velocity by making (2/(0) = fin- Note that here a prime 
indicates derivation with respect to the compactified coordinate x. 

At x = 1 we impose asymptotic flatness by requiring 

m(l) = 1 
n( 1) = 1 

/(l) = 1 (62) 

w(l) = 0 
a v ( 1 ) = 0 . 

Alternatively, we can also impose the angular momentum by requiring u>'( 1) = — J/(47t 2 ). 

In the system of equations we are using, the a 0 function was eliminated by introducing a first integral of the system 
in terms of the electric charge. Using these boundary conditions we can obtain configurations with selected pairs 
of values for (Q , Hr) or (Q, J) for a given value of A. Once a configuration is obtained, we can compute the global 
charges. 

From the asymptotic behaviour of the functions at x = 1 we extract the mass M, the angular momentum J, and 
the magnetic moment g m ag- With these quantities we calculate the gyromagnetic factor g. From the behaviour of 
the functions at the horizon i = 0 we extract the horizon mass Mu, the horizon angular momentum Jh, the horizon 
area Ah and the horizon electrostatic potential <£■#. The Srnarr formula is always satisfied within the precision of the 
black hole configurations computed. 


V. GLOBAL SOLUTIONS: NUMERICAL RESULTS 

We now discuss the global EMCS black hole solutions obtained by numerical integration as discribed above, with 
emphasis on extremal configurations. Here we focus on the most interesting case for the CS coupling constant A, 
namely A > 2. We first discuss the relation between the near-horizon solutions and the global solutions, and exhibit 
the intriguing branch structure of the latter. Subsequently we address the sequence of J = 0 solutions and reveal the 
node structure of the two functions a v {r) and w(r). Finally, we discuss the domain of existence of the global solutions. 

Concerning the lower values of A, let us remark that for 0 < A < 1.91 the global extremal solutions are in one to 
one correspondence to the near-horizon solutions, while in the transitions region 1.91 < A < 2 this does no longer 
seem to be the case. When approaching the RN solution for values of A in the transition region and positive values of 
the charge, the functions of the rotating global extremal solutions develop highly divergent derivatives. This forbids 
definite statements on the existence and properties of the numerically constructed global solutions in this area of 
parameter space. However, our analysis indicates, that the rotating global extremal solutions are disconnected from 
the static extremal RN black hole in the transition region. In particular, we did not obtain any non-static extremal 
J = 0 solutions for A < 2. 


A. Near-horizon versus global solutions: the branch structure of extremal black holes 

Let us now recall the near-horizon solutions and overlay the corresponding properties of the global solutions for 
direct comparison. In Fig. [4] we exhibit the horizon area Ah and the horizon angular momentum Jh versus the total 
angular momentum J for both sets of solutions. For definiteness, we have chosen CS coupling Acs = 5 and charge 
Q = ±1. However, the pattern observed has the same structure for all values of the CS coupling constant A > 2 and 
values of the charge Q. The near-horizon solutions are marked by thick lines, while the global solutions are marked 
by thin lines in the figures. 

Clearly, for the negative charge solutions, shown in red, both global and near-horizon solutions match. Thus 
in this case all near-horizon solutions are realized globally. However, for positive charge solutions, shown in green 
(near-horizon) and blue (global), the situation is very different. The branches of global solutions end in the cusps 
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FIG. 4: Global versus near-horizon solutions: The horizon area Ah (a) and the horizon angular momentum Jh versus the 
angular momentum J for charge \Q\ = 1 and for CS coupling A = 5. The thick lines represent the near horizon solutions 
for positive charge (green) and negative charge (red), while the thin lines the global solutions for positive charge (blue) and 
negative charge (red). The cusps are marked by x and denoted by Cr i, the bifurcation points are marked by + and denoted 
by -£>[...] (see text). The black asterisk represents the extremal RN solution. 


C[i j2 ], while the near-horizon solution branches extend further to C n h- Thus the near-horizon solutions between both 
cusps are not realized globally. Moreover, the full near-horizon branch connecting the cusp C n h and the extremal RN 
solution, which is marked by an asterisk, is neither realized globally - with the exception of the RN solution itself, of 
course. 

Therefore we conclude, that the static extremal RN solution with Q > 0 is rotationally isolated from the family of 
global stationary extremal solutions, when A > 2 (and presumably already for A > 1.91). The horizon area Ar of the 
global solutions with J = 0 is significantly smaller than the horizon area of the RN solution. Moreover, they carry 
finite horizon angular momentum Jh. 

We note, that the existence of near-horizon solutions, which do not possess global counterparts, has been noted first 
by Chen et al. [36] for the extremal dyonic black holes of D = 4 Gaufi-Bonnet gravity. However, in the EMCS case 
the relation between near-horizon solutions and global solutions is even more surprising, as discussed in the following. 
By considering the physical quantities, not accessible in the near-horizon formalism, like the mass M or the horizon 
angular velocity IIh, we now present the intricate pattern of branches of the global rotating extremal black holes. 

We exhibit in Fig. [5] the mass M (a) and the horizon angular velocity fin (b) versus the total angular momentum J 
of the global extremal black hole solutions. We have chosen CS coupling constant A = 5 and electric charge Q = ± 1. 
These figures contrast the simple branch structure for negative charge (red) with the highly complicated branch 
structure for positive charge (blue). _ 

Let us begin the detaile d dis cussion of the branch structure for po sitive Q with Fig. |5(a)[ which has a number 
of similar features to Fig. |4(a)| For instance, the two points in Fig. |5(a)| where a minimal value of the mass is 
encountered, correspond to the two singular solutions in Fig. |4(a)| which have vanishing horizon area. We now choose 
one of these particular solutions, the one with J < 0, as the starting point of our analysis of the branch structure. 

When we decrease the angular momentum from this Jh < 0, An = 0 solution, we obtain a single infinite branch of 
extremal solutions, whose mass increases monotonically with decreasing J. In contrast, when we increase the angular 
momentum from this Jh < 0, An = 0 solution, we cross the first non-static J = 0 solution, which we label n = 1. 
Increasing J further, we then cross a particular solution, which we label JJ[i 2 ] - Eventually, this branch ends at a finite 
maximal value of J, where a cusp is encountered, which we label C n 2 i ■ The solution at the cusp is perfectly regular. 

From the cusp Co , 2 ] a new branch of solutions bends backwards towards smaller masses, as the angular momentum 
is decreased. Along this branch first another particular solution is encountered, which we label _B[ 2 , 3 ]- Then the 
branch crosses the next non-static J = 0 solution, labeled n = 2. It eventually ends in a solution labeled 2 j, when 
it reaches the J —> — J symmetric first branch at the solution -B[i, 2 ]- Thus, we have encountered a bifurcation point, 
where the two branches meet, which has two distinct solutions, /?,* 2 , and that possess the same global charges. 
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(a) (b) 


FIG. 5: Global solutions: The total mass M (a) and the horizon angular velocity S!h (b) versus the angular momentum J 
for charge Q = 1 (blue), Q = —1 (red) and for CS coupling A = 5. The cusps are marked by x and denoted by Cr i, the 
bifurcation points are marked by + and denoted by By i (see text). The non-static J = 0 solutions are numbered by n and 
marked by red dots in different colors. The black asterisk represents the extremal RN solution. 



J 


FIG. 6: Schematic representation of the branch structure of the global solutions: The total mass M versus the angular 
momentum J. The cusps are marked by x and denoted by Cri, the bifurcation points are marked by + and denoted by By i. 
The non-static J = 0 solutions are numbered by n. 


This reveals the presence of non-uniqueness among the extremal solutions. While the global charges of these 
solutions are the same, and their horizon angular velocities are also the same, as seen in Fig. |5(b)| they differ in their 
horizon area, as seen in Fig. 4(a) The area of B* ± 2 j is considerably smaller than the area of By 12 ]- 
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These two sets of symmetric branches with cusps C[i, 2 ] and bifurcation points B f\ 2 j form only the lowest 

mass solutions of a whole tower of branches, built in an analogous ever repeating manner. We have already noted the 
presence of the particular (symmetric) points B[ 2 ,3]- Here the next (symmetric) bifurcations are encountered, where 
new branches arise, whose initial solutions are labeled B* 2 3 j. 

Starting at the solution B* 2 3 j with negative J , we cross the next non-static J = 0 solution, n = 3. Increasing J 
further, we then cross the solution at the next bifurcation point which we label B[ 3 M, and then encounter the next 
cusp C[ 3 ) 4 ]. From CA 41 a new branch of solutions bends backwards towards smaller masses, as the angular momentum 
is decreased. Now the next bifurcation point is encountered, _B [4 51 , then the next non-static J = 0 solution, n = 4, is 
passed, and the branch ends at the bifurcation point, formed by B * 3 4 , and H [ 3 4 ]. 

This general scheme is exhibited schematically in Fig. [ 6 ] Starting from the solution B^ k+1 fc+2 j with negative J, 
first the (k + 2)th non-static J = 0 solution is encountered, then the bifurcation point B[ k+2 ^ k + 3 ], and subsequently 
the cusp C[fe + 2 i fc+ 3 ] • From C\j~+ 2 ,k+ 3 \ a new branch bends backwards towards smaller masses, encountering first the 
bifurcation point B^ k+3 k+4 ] then the (k + 3)th non-static J = 0 solution, and it ends at the bifurcation point, formed 

b y B [k+2,k+3] and B [k+2,k+3]- 

For any fc, the two solutions at the bifurcation points B[ k+2k+3 ] and B * k+2 fc+3] possess the same global charges, 
whereas their horizon area and horizon angular momentum differ. Thus these solutions constitute a sequence of 
extremal solutions, that violate uniqueness. The horizon angular velocities of the solutions B^ k+2 ,k+ 3 ] and B* k+2 fe+3 |, 
however, possess the same value. This is clearly seen for the lowest k in Fig. |5(b)| This figure also highlights the 
singular nature of the points with vanishing area, showing that the horizon angular velocity jumps from to — fin 
at these points. 

Having generated in this way a whole sequence of branches, labeled by n, let us inspect the extent of these branches. 
Here we observe that with increasing n, the cusps C njn + 1 occur at decreasing values of | J\. Thus the range of extremal 
global solutions for fixed | J| decreases with increasing n. 

Coming finally back to the comparison of global and near-horizon solutions, we conclude that a given near-horizon 
solution can correspond to 

i) more than one global solution , 

ii) precisely one global solution , 

iii) no global solution at all. 

In fact, we conjecture, that a given near-horizon solution may even correspond to an infinite set of global solutions, 
as discussed in the following. 


B. Radial excitations: J = 0 solutions 


Let us now address the question of how these branches of global extremal solutions differ. In particular, we would 
like to associate the integer n counting the branches with a physical property of the solutions along these branches. 
To that end, let us first focus on the set of J = 0 solutions. As seen in Fig. |5(a)| the mass M n of each pair of 
degenerate J = 0 solutions increases with n. This suggests, that the solutions form a sequence of higher and higher 
excited states, with the lowest mass solution forming the ground state or fundamental solution. 

Excited states of otherwise very similar solutions are typically found, when radial excitations are allowed. Examples 
range from the hydrogen atom, via sphaleron solutions in flat space (37}, to hairy black holes 381, all featuring an 
infinite number of radial excitations of the ground state. 

Indeed, it is the occurrence of radial excitations, which results in the sequence of J = 0 solutions and the associated 
intriguing pattern of branches of the global extremal black holes. We exhibit in Fig. [7] the gauge field function a^ 
(a) and the metric function w (b) versus the radial coordinate. Again we have chosen charge Q = 1 and CS coupling 
constant A = 5. 


As the radial coordinate in the figure we have employed the logarithm of the scaled radial coordinate log 10 . 

Recall, that all these solutions possess the same horizon area Ah, as seen in Fig. |4(a)| The first node always refers 
to spatial infinity. We note, that in this logarithmic scale, the nodes of the solutions are located at roughly equal 
spacings. Moreover, the kth node of the excited solutions with more than k nodes are located very close to each other. 
(This is also known for the hydrogen atom and other radial excitations.) 

For the metric function w we have chosen to exhibit the logarithm of its modulus, since in contrast to a v it varies 
over many orders of magnitude, as n is increased. However, as seen in Fig. |7(b)| beyond the nodes, for small values 
of the radial coordinate, the functions log 10 |o;| exhibit the same slope for all n. This also holds for large values of the 
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l°gl0( r ^H ) 
(a) 



FIG. 7: Global solutions: The gauge field function a v (a) and the logarithm of the modulus of the metric function to (b) versus 
the logarithm of the scaled radial coordinate log 10 ^ r/A . Shown are the lowest radial excitations n = 1,7 of a v (a) and 
n = 1,10, 15 and 31 of to (b) for charge Q = 1 and CS coupling constant A = 5. 


radial coordinate. Besides the ten lowest n, Fig. |7(b)| also exhibits the metric function for n = 15 and n = 31. This 
clearly shows, that the pattern continues to high values of n. 

Solutions with n in the thirties are the solutions with the highest node numbers we could achieve numerically so far. 
However, we conjecture, that this sequence of J = 0 solutions can be continued to arbitrarily high values of n. Thus 
we conjecture, that we have obtained the lowest members of an infinite sequence of radially excited J = 0 solutions. If 
this is really the case, then the corresponding near-horizon solutions will correspond to an infinite sequence of global 
solutions. 

As mentioned already, with increasing node number n, the mass M„ increases monotonically. In fact, it converges 
monotonically from below to the mass of the extremal static RN black hole, Mrn- The Smarr formula then implies, 
that the horizon electrostatic potential (fin converges likewise to the extremal static RN value, ■Frn- At the same time, 
the horizon angular velocity Q n converges to zero, and thus also to the extremal static RN value. The horizon area 
Ah, however, is independent of n and has the same value for all the members of the (presumably) infinite sequence, 
and differs significantly from the extremal static RN value. The same holds true for the magnitude of the horizon 
angular momentum Jh- 

Let us address now the A-dependence of the observed pattern of extremal J = 0 solutions. As pointed out already, 
there is a lower critical value A cr , beyond which no rotating J = 0 solutions exist. When this critical value is 
approached from above, the mass of the rotating J = 0 solutions should tend to the mass of the extremal static RN 
solution, in order to have a smooth emergence of these solutions at the critical value. 

In Fig. [8] we demonstrate that this is indeed the case. Here we present the mass M of the rotating J = 0 solutions 
with charge Q = 1 and up to seven nodes versus the CS coupling A. As expected the mass M n of the set of solutions 
with n nodes approaches the mass of the static extremal RN solution Mrn, when the CS coupling tends to A = 2. In 
particular, we observe that the larger the value of n, the closer is the mass M n to the limiting mass Mrn for fixed 
A. To discern the branch structure, it is therefore preferable to employ values of A not too close to the critical value. 
This explains our choice A = 5 in most figures. 

In Fig. Owe demonstrate the sequence of J = 0 solutions for a higher value of the CS coupling constant, A = 10. 
Obviously, the figures for the gauge function and the metric function are completely analogous to their counterparts, 
exhibited in Fig. [7| 

The nodal excitations seen in the sequence of J = 0 solutions, are also present in the general set of global extremal 
solutions. Along the various branches of solutions, new nodes appear from or disappear towards radial infinity. Thus 
all global extremal solution possess an associated node number. However, the extent of the branches decreases with 
increasing n. Thus for a given J / 0, only a finite number of radial excitations may exist. 

Likewise, also the non-extremal solutions possess an associated node number. Fixing, for instance, the horizon 
radius, families of non-extremal solutions are obtained, which exhibit an intricate branch structure reminiscent of the 
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FIG. 8: Global solutions: The total mass M versus the CS coupling constant A. Shown are the lowest radial excitations 
n = l,7 for charge Q = 1. 



l°gio<hA H ) 



(a) 


(b) 


FIG. 9: Global solutions: The gauge field function a v (a) and the logarithm of the modulus of the metric function to (b) versus 
the logarithm of the scaled radial coordinate log 10 ^ r/A . Shown are the lowest radial excitations n = 1,..., 7 of a v (a) and 
n = l,..., 10, 15 and 32 of to (b) for charge Q = 1 and CS coupling constant A = 10. 


one of the extremal solutions, but with less branches. However, the number of branches increases the more, the closer 
the family of non-extremal solutions approaches the set of extremal solutions. 

Finally let us remark that the node number might possibly be considered in the quest to regain uniqueness. 


C. Domain of existence for A > 2 

Finally, we would like to address the domain of existence of global EMCS black hole solutions. Here we would 
like to first consider the global charges and, in particular, the angular momentum versus the charge for fixed mass. 
Typically, the boundary of the domain of existence of black hole solutions is then formed by extremal solutions, 
where the non-extremal solutions reside inside this boundary. Therefore we exhibit in Fig. |10(a)| the scaled angular 
momentum j = J/M 3 / 2 versus the scaled charge q = Q/M for extremal solutions, choosing again CS coupling A = 5. 
(A similar picture is found for other values of A > 2.) The figure includes the first few bifurcation points, cusps and 
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FIG. 10: Global and horizon properties: The scaled angular momentum j = J/M 3 ^ 2 of extremal solutions including the 
bifurcation points and cusps versus the scaled charge q = \Q\/M (a). The scaled angular momentum j = J/M 3//2 of extremal 
and selected sets of non-extremal solutions (b). The scaled horizon angular velocity Oh M 1 ^ 2 (c), and the scaled surface gravity 
kM 1 / 2 (d) for the same sets of extremal and non-extremal solutions. (CS coupling A = 5) 


non-static J = 0 solutions. This allows us to connect the locations of the solutions in the domain of existence with 
their locations in the previous figures. 

Fig. |10(a)| then shows, that when the scaled angular momentum is considered versus the scaled charge, the domain 
of existence is indeed delimited by branches of extremal solutions, which constitute its boundary. Clearly, the domain 
of existence is different for positive and negative charge solutions, since for spinning solutions the CS term breaks the 
charge reversal symmetry, Q —> — Q. The extremal solutions with negative charge form the smooth Q < 0 boundary. 
This contains the extremal static RN solution, when the angular momentum vanishes. The extremal solutions with 
positive charge, on the other hand, reflect the complicated branch structure discussed above. 

For positive charge only a subset of the extremal solutions forms the Q > 0 boundary. This subset of solutions 
starts at the extremal Myers-Perry (q = 0) solutions, passes the cusps with vanishing horizon area, and ends at the 
first (n = 1) non-static J = 0 solution. Thus in contrast to the Q < 0 boundary, the q > 0 boundary is not smooth, 
since it contains the pair of symmetric J —> — J cusps, which reside at the maximal value q mayi of the scaled charge q. 
The solutions at the cusps correspond to the singular solutions with vanishing horizon area in Fig. 0] Moreover, this 
maximal value of q is larger than the value of q of the corresponding extremal RN solution. 

Even more surprising, however, is the fact, that there are extremal solutions in the interior of the domain of 
existence. In particular, the Q > 0 extremal RN solution is no longer part of the boundary. Such a feature does not 
occur in theories like pure Einstein-Maxwell theory, or Einstein-Maxwell-dilaton theory, and neither for EMCS theory 
for sufficiently low values of the CS coupling constant. This intriguing fact produces a peculiar type of uniqueness 
violation: there exist both extremal and non-extremal solutions, which possess the same global charges, M, J, and 

Q. 
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For more clarity and later reference, we replicate the demonstration of the domain of existence in Fig. |10(b)| now 
without the inset, the bifurcations and the cusps, but instead with the inclusion of a number of selected sets of non¬ 
extremal solutions. These sets contain, in particular, the branches of non-static solutions with non-rotating horizon, 
= 0. Clearly, all solutions lie within the domain of existence, while the crossings of lines exhibit the non-uniqueness 
of the solutions. 

Let us next consider the horizon properties of the solutions. In Fig. 10(c) we exhibit the scaled horizon angular 
velocity fln-M 1 / 2 for the extremal solutions and the same set of non-extremal solutions. Most of the boundary for 
the scaled horizon angular velocity fln-M 1 / 2 is formed by a subset of the extremal black holes, only the righthand 
boundary, i.e., the large q limit, is given by the maximum value of q , q max . All other extremal and all non-extremal 
solutions reside within these boundaries. 

We exhibit the scaled surface gravity kM 1 / 2 in Fig. |10(d)| Clearly, the extremal black holes form the lower part 
of the boundary, since they all possess re = 0. For the Q < 0 solutions, the static black holes form the upper part of 
the boundary, as is often the case. However, for the Q < 0 solutions, the static solutions represent only the first part 
of the upper part of the boundary. Then another set of solutions adopts the role of forming the upper part of the 
boundary. Interestingly, this set corresponds precisely to the first set of non-static Hh = 0 solutions, which start at a 
certain static RN solution and reach all the way to the singular cusp. 

Let us next address the scaled area an = Ah/M 3 / 2 . For Q < 0, the extremal and the static solutions form again 
together the boundary. For Q > 0, the upper boundary part is also formed the first part of the set of static solutions, 
until the first set of non-static fin = 0 solutions emerges from the static solutions, which then form the remaining 
upper part of the boundary. 

Finally, we address the gyromagnetic ratio g of the same sets of black hole solutions. Again, for Q < 0, the 
gyromagnetic ratio is bounded by the static and extremal solutions. For Q > 0, however, the gyromagnetic ratio 
becomes unbounded, as was noted before 


VI. CONCLUSIONS 


Here we have discussed global and near-horizon solutions of EMCS theory in five dimensions, focussing on solu¬ 
tions with equal magnitude angular momenta, to enhance the symmetry of the solutions, making the analytical and 
numerical analysis much more tractable while at the same time revealing already numerous intriguing features of the 
solutions. 

Since the CS term breaks the charge reversal invariance, the families of solutions for negative charge and positive 
charge no longer agree, when the CS term contributes, be., for rotating solutions. Choosing the CS coupling constant 
A to be positive, then the negative charge solutions do not represent any particular peculiarities. The positive charge 
solutions, however, develop a number of very interesting properties, as the CS coupling constant A first reaches the 
supergravity value Asg and then exceeds it. 

At Asg the solutions are known analytically [22j, with the BMPV solutions [20!] forming a special case, which 
represent extremal non-static Oh = 0 solutions, ending in a singular Ah = 0 solution. For A / Asg the global 
solutions are known only perturbatively EMI or numerically |23l425| . obtained by solving the corresponding set 
of ordinary differential equations, subject to a set of appropriate boundary conditions. However, the near-horizon 
solutions can be obtained in a simpler way, by making use of an appropriately modified near-horizon formalism, where 
the modification is caused by the CS term. 

The presence of a CS term always implies the occurrence of (a set of two degenerate J —> — J symmetric) singular 
Ah = 0 solutions. But only for A > Asg these singular solutions appear prominent in the domain of existence, residing 
at a cusp. We exhibit the A-dependence of these solutions in Fig. [TlJ The cusp points have been obtained both for 
the near-horizon and the global solutions, as demonstrated in Fig. 11(a) When the cusp points are considered for the 
scaled angular momentum and the scaled charge, i.e., in the form j(q), we note, that they form an almost straight 
line, as depicted in Fig. |ll(b)| 

Here we have focussed our considerations on values of the CS coupling constant A > 2 Asg> since above this limit a 
new type of extremal solutions is present: extremal non-static J = 0 solutions. These solutions come as degenerate 
pairs, symmetric with respect to J —> — J, where the lowest solution resides in the boundary of existence of the black 
hole solutions. The solutions form a presumably infinite sequence, that can be labeled by an integer n. This integer 
is associated with the number of nodes of the metric function oj and the gauge field function a v . To our knowledge, 
this is the first example of black holes with Abelian fields which form excited states, which are reminiscent of the 
radial excitations of atoms, of sphalerons, or of hairy black holes. 

With increasing node number the mass of these extremal non-static J = 0 black holes converges to the mass of the 
extremal static RN black hole, while their horizon angular velocity tends to zero along with their magnetic moment. 
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(a) 


(b) 


FIG. 11: Singular Ah = 0 solutions: The scaled angular momentum J/Q 3 ^ 2 versus the CS coupling constant for near-horizon 
and global solutions (a); the scaled angular momentum j = J/M 3 ^ 2 versus the scaled charge q = \Q\/M for global solutions 
(between the MP solution and the Ah = 0 cusp) and several values of A. 


Their area, however, retains the same constant value, which is different from the static value. We demonstrate this 
pointwise convergence of the solutions in the Appendix. 

The non-static J = 0 solutions are located symmetrically within an intriguing pattern of branches of extremal 
solutions. Most of these branches of extremal solutions reside within the domain of existence of EMCS black holes, 
when the scaled angular momentum is considered versus the scaled charged. This is rather unusual, since the extremal 
solutions typically form the boundary of the domain of existence in this case. It has the interesting consequence, 
that there is non-uniqueness between extremal and non-extremal black holes. Furthermore, there is non-uniqueness 
between extremal black holes, which reside at the bifurcation points of the branches, while non-uniqueness between 
non-extremal black holes was reported already earlier [23|. 

The comparison between near-horizon solutions and global solutions has also led to surprises, namely we have seen, 
that a given near horizon solution can correspond to i) more than one global solution , possibly even an infinite set , 
ii) precisely one global solution, or iii) no global solution at all. Thus our findings show that the intuition based on 
known exact solutions cannot be safely applied in the general case. Moreover, the results of a near-horizon analysis 
of extremal solutions cannot be always safely extrapolated to the global case, as assumed sometimes in the literature. 

Similar black holes can be obtained in EMCS theory in higher odd dimensions, when all the angular momenta have 
the same magnitude (24| . In particular, our preliminary results show that for D = 7,9 this sequence of non-static 
J = 0 solutions is also found, and the solutions can be characterized by the number of nodes of the corresponding 
metric and gauge field functions. 

We conjecture that extremal black holes with simil ar p roperties may also exist in other theories, in particular, in an 
Einstein-Maxwell-dilaton theory in four dimensions [39j |. In this case, however, a set of partial differential equations 
must be solved, making the corresponding analysis much more involved. 

However, our next step will be the inclusion of a cosmological constant for these EMCS black holes. When the new 
intriguing phenomena observed for the asymptotically flat solutions survive in the presence of a negative cosmological 
constant, the AdS/CFT correspondence may imply interesting consequences for the associated four-dimensional field 
theories, living on the boundary. 
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VII. APPENDIX 

In this appendix we briefly illustrate the full set of functions and the dependence of the nodes for the sequence of 
non-static J = 0 extremal solutions and compare with the corresponding RN functions. 
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FIG. 12: Global J = 0 solutions: The metric function / (a), m (b), n (c), uj (d) and the gauge field functions ao (e) and a v (f) 
versus the radial coordinate for the lowest radial excitations n = 1,7. For comparison, the corresponding functions for the 
extremal RN solution are also shown. (Charge <5 = 1 and CS coupling A = 5). 


We present in Fig. [T2] the radial dependence of the functions of the fundamental J = 0 solution and the first few 
excitations. We note that the functions /, m, n and ao exhibit a fast convergence towards a limiting solution, which 
is represented by the static extremal RN solutions. However, this convergence is only local. At the horizon, x = 0, 
the functions do not tend to the RN values in most cases. In the limit n —> oo the functions will then exhibit a 
jump. We notice this also for the function a v . Such a discontinuity is necessary, to guarantee that these solutions, 
while approaching the corresponding RN solution, retain a horizon area which differs from the RN horizon area, but 
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remains the same for all solutions of the sequence. 
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